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SUMMARY 


Two  methods  of  using  the  concept  of  linearized  characteristics  are 
derived  for  the  one-dimensional  unsteady  flow  in  a  tube  that  is  rotated 
about  an  axis  perpendicular  to  the  axis  of  the  tube.  One  of  the  methods 
corresponds  to  that  used  by  Ferri  in  his  basic  work  on  the  subject. 
Solutions  are  made  by  both  methods  for  boundary  conditions  that  allow 
analytic  solutions.  Comparison  shows  that  both  methods  give  the  same 
results  but  there  are  significant  differences  in  their  application. 


INTRODUCTION 


The  term  linearized  characteristics  has  been  applied  by  Ferri  in 
reference  1  to  a  process  of  superposing  a  small  perturbation  on  a  basic 
solution  of  a  set  of  nonlinear  hyperbolic  differential  equations.  The 
perturbation  can  be  due  to  a  change  in  the  prescribed  conditions  along 
the  boundary ,  in  the  position  of  the  boundary,  or  in  the  differential 
equations  themselves.  The  unknowns  can  each  be  expanded,  for  example, 
into  a  power  series  in  a  parameter  representative  of  the  perturbation, 
with  unknown  coefficients,  except  for  the  coefficient  of  the  first  term, 
representing  the  basic  solution,  which  is  known.  The  equations  for  the 
other  coefficients  are  found  to  be  linear  and  can  be  solved  by  charac¬ 
teristic  methods.  As  applied  to  steady  supersonic  flow,  this  procedure 
has  the  advantage,  stated  roughly,  over  ordinary  perturbation  theory 
that  it  can  be  used  for  perturbations  of  any  known  flow  (to  which  the 
characteristic  method  of  solution  can  be  applied)  rather  than  only  for 
perturbations  of  constant  parallel  flow. 

A  version  of  the  process  was  used  by  Sauer  (ref.  2)  for  determining 
the  steady  supersonic  flow  about  bodies  of  revolution  at  small  angles  of 
attack  but  apparently  has  not  been  used  for  other  problems  except  in 
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reference  1 ,  where  its  application  to  a  variety  of  problems  is  discussed. 
In  references  1  and  2,  the  expansion  of  the  unknowns  was  substituted 
into  the  original  equations  of  motion,  which  were  then  transformed  into 
characteristic  form.  The  purpose  of  the  present  paper  is  to  compare 
that  method  of  procedure  (called  method  A  subsequently)  with  an  alter¬ 
nate  method  (method  B)  whereby  expansions  of  the  dependent  and  inde¬ 
pendent  variables  are  substituted  into  the  characteristic  form  of  the 
equations . 

The  comparison  is  made  by  developing  the  equations  of  both  methods 
as  they  apply  to  a  simple  set  of  nonlinear  hyperbolic  equations  of  gas 
dynamics,  the  set  relating  to  the  one -dimensional  unsteady  motion  of  a 
gas  in  a  tube  of  constant  cross  section  with  the  addition  of  a  perturba¬ 
tion  due  to  the  rotation  of  the  tube.  The  equations  derived  for  this 
case  are  then  solved  for  a  set  of  boundary  conditions  for  which  the 
solution  can  be  found  analytically  by  both  methods. 


a, A 
cv 

k 

P 

R 

s, S 

t, T 

u, U 
x,X 
<x,P 
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SYMBOLS 

speed  of  sound 

specific  heat  at  constant  volume 

unknown  functions  of  one  variable 

constant 

pressure 

gas  constant 

ehtropy  increment  divided  by  gas  constant 
time 

velocity 

distance  along  tube  measured  from  axis  of  rotation 
characteristic  parameters 
ratio  of  specific  heats 
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absolute  temperature 
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R  7-1 

P  density 

angular  velocity 

Subscripts: 

Letter  subscripts  indicate  partial  derivatives. 

Number  subscripts  indicate  order  of  approximation  with  the  following 
exceptions : 

a0  reference  value  of  a 

Xq  and  xp  inner  and  outer  ends  of  tube,  respectively 

ANALYSIS  OF  UNSTEADY  ROTATING  FLOW 
Development  of  Equations 


The  linearized-characteristics  concept  is  applied  to  a  relatively 
simple  problem  in  order  to  bring  out  more  clearly  some  of  its  essential 
features.  Consider  the  unsteady  motion  of  an  inviscid  gas  in  a  cylin¬ 
drical  tube  which  can  rotate  with  angular  velocity  co  about  an  axis 
perpendicular  to  the  generators  of  the  cylinder  (fig.  1).  Distance 
along  the  tube  is  denoted  by  x  measured  from  the  axis  of  rotation. 
Only  values  of  x  large  in  comparison  with  the  greatest  lateral  dimen¬ 
sion  of  the  tube  in  a  plane  of  rotation  are  considered  so  that  the 
local  velocity  of  the  rotation  may  be  taken  as  cox.  Cross  flows  and 
lateral  pressure  gradients  are  neglected  so  that  the  flow  considered 
is  one-dimensional.  A  perfect  gas  with  constant  specific  heats  is 
assumed. 


When  co=0,  the  ordinary  equations  of  motion  apply  (ref.  3):  that 
is,  the  continuity  equation 


pt  +  Pux  +  upx  =  0  (!) 


ut  +  uux  -  "  P  Px  (2) 


and  the  momentum  equation 
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When  a)  ^  0,  equation  (l)  is  unaltered,  but  equation  (2)  has  another 
term  representing  the  acceleration  along  the  tube  due  to  the  rotation 

ut  +  uux  -  co2x  =  -  Px  (3) 

The  entropy  is  now  assumed  to  be  constant  throughout  the  gas  so  that 

P  =  kp7  (k) 

See  appendix  A  for  a  treatment  of  the  case  of  variable  entropy.  Equa¬ 
tion  (k),  the  equation  of  state 

p  =  pR0 


and  the  relation 


a2  =  ^  =  7R0 
dp 


allow  the  replacement  of  p  and  p  by  functions  of  a  in  equations  (l) 
and  (3).  Thus  the  equations  of  motion  become 

Xa^  +  aux  +  Xuax  =  0 

)>  (5) 

p 

u^.  +  uux  +  Xaax  -  co  x  =  0 

where 

x  =  =  2 

^  R  7-1 

These  equations  are  referred  to  subsequently  as  the  "primary  equations 
to  distinguish  them  from  the  characteristic  equations  which  are  now 
developed. 

Adding  and  subtracting  equations  (5)  gives 

Uj.  +  (u  +  a)^  +  xjat  +  (u  +  a)a^j  -  co2x  =  0 

f  (6) 

ut  +  (u  -  a)ux  -  xjat  +  (u  -  a)a^j  -  co2x  =  0 

By  the  usual  procedure  for  hyperbolic  equations  (ref.  3)>  new  independent 
variables  a  and  (B  are  introduced  so  that  one  equation  of  equations  (6) 
has  derivatives  with  respect  to  a  only  and  the  other,  with  respect 
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to  P  only.  Then  the  lines  a  =  constant  and  p  =  constant  form  the 
two  sets  of  characteristics.  This  process  is  accomplished  by  writing 


so  that 


xa  =  (u  +  a)ta 
Xp  =  (u  -  a)tp 


s_ 

5a 


(u  +  a)~ 

OX 


(7) 


5  _  + 

—  =  tp 

5p 


—  +  (u  -  a)  — 
5t  5x 


Thus,  equations  (6)  become 


Uq,  +  Xaa  -  a/xta  =  0 


Up 


-  \aR  -  aftctR  =  0 


(8) 


Equations  (7)  and-  (8)  form  a  set  of  four  equations  in  the  four  unknowns 
u,  a,  x,  and  t  and  are  equivalent  to  equations  (5)*  These  equations 
can  always  be  solved  by  step-by-step  methods. 


Method  A.-  Linearized  approximations  to  equations  (5)  are  now 
developed.  Let 

.  ^ 
u(x,t)  =  U  +  ClAlj  +  CJO^Up  +  .  .  . 

> 

a(x,t)  =  A  +  co2a1  +  +  .  .  . 


(9) 


Only  even  powers  of  o>  are  used  since  equations  (5)  show  that  the  problem 
is  independent  of  the  sign  of  <jq.  By  substituting  expressions  (9)  into 
equations  (5)  and  setting  the  coefficients  of  each  power  of  00  equal  to 
zero,  the  following  equations  are  obtained: 


XAt  +  AUX  +  XUAx  =  0 

l 

Ut  +  UUX  +  XAAx  =  0 


(10) 
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+  Xa2a  +  (Ux  +  XAx)ta(u2  +  a2)  +  (ulx  +  Aalx)ta(ul  +  al)  =  0 

>  (16) 

ugp  -  Xa2p  +  (Ux  -  XAx)tp(u2  -  a2)  +  (ulx  -  Xalx)tp(u1  -  aX)  =  0 
and  so  forth. 

The  development  and  use  of  equations  (13)  to  (l6)  is  referred  to 
hereafter  as  "method  A." 

Method  B . -  Instead  of  substituting  expansions  of  u  and  a  into 
the  primary  equations  (eqs.  (5))  as  in  method  A,  these  expansions  are 
now  substituted  into  the  characteristic  form  of  the  equations  (eqs.  (8)). 
But  the  values  u  and  a  must  be  considered  functions  of  a  and  3 
in  equations  (8),  and  the  corresponding  values  of  x  and  t  are  found 
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from  equations  (7).  Hence,  x  and  t  are  also  expanded.  The  method 
is,  therefore,  similar  to  the  technique  presented  by  Lighthill  in  refer¬ 
ences  4  and  5* 

By  writing 

u(a,|3)  =  U  +  a>2Up  +  oAj2  +  .  .  . 
a(a,|3)  =  A  +  a>2Ap  +  m^A2  +  .  .  . 


x(a, 3)  =  X  +  ay^Xp  +  0^X2  +  .  .  . 

t(a, 3)  =  T  +  cn2Tp  +  +  .  .  . 

substituting  into  equations  (7)  and  (8),  and  setting  the  coefficient  of 
each  power  of  cn  equal  to  zero,  the  following  sets  of  equations  are 
obtained: 


I 

; 


(17a) 


(17b) 


Ua  +  -  0 


Up  -  XAp  =  0 


(18a) 


Xa  -  (U  +  A)Ta  =  0 
Xp  -  (U  -  A)Tp  =  0 


Ula  +  ^la  = 
U13  "  ^lp  = 


(18b) 


(19a) 


Xia  -  (U  +  A)Tla  =  (Up  +  Ap)Ta 
Xpp  -  (U  -  A)Tpp  =  (Up  -  Ap)Tp 


(19b) 


u2a  +  XA2a 
U2p  -  AA2p 


=  XlTa  +  XTla 


=  X1TP  +  XT1P 


(20a) 
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X2ot  -  (U  +  A)T2a  =  (Ui  +  Ai)Tia  +  (U2  +  A2)Tal 


X2p  -  (U  -  A)T2f3  =  (Ui  -  A!)Tip  +  (U2 


(20b) 


and  so  forth.  The  development  and  use  of  these  equations  is  referred 
to  as  "method  B." 


Equations  (18)  and  equations  (13)  and  (l4)  are  exactly  the  same 
and  their  solution  is  the  complete  solution  of  the  problem  when  cn  =  0. 
This  solution,  the  basic  solution,  is  essential  to  the  solution  of  the 
higher  approximations  and  is  considered  known.  The  equations  -for  the 
higher  approximations  in  both  methods  are  all  linear  in  the  unknowns. 


1 

Application  to  Specific  Boundary  Conditions 

In  order  to  demonstrate  some  of  the  differences  between  the  two 
methods,  an  application  to  a  simple  case  for  which  the  equations  can 
be  integrated  analytically  is  now  made.  Assume  that  the  tube  extends 
from  x  =  x0  to  x  =  x^_  (X2  >  xq)  and  is  closed  at  the  ends.  The 

tube  is  initially  at  rest,  and  the  state  of  the  gas  is  given  by  u  =  0 
and  a  =  aQ,  a  constant.  At  time  t  =  0  the  tube  is  started  rotating- 

about  x  =  0  with  angular  velocity  o>.  The  neglect  of  cross  flows  and 
lateral  pressure  gradients  might  be  completely  unjustified  in  this  case; 
nevertheless,  the  problem  is  interesting  mathematically. 

The  problem,  then,  is  to  solve  the  primary  equations  (eqs.  (5)) 
with  the  following  boundary  conditions: 

At  t  =  0,  ' 

u  =  0;  ’  a  =  aQ  (x0  x  <  xq)  (21a) 


At  x  =  x0  and  x  =  xq_. 


u  =  0  (t  >  0) 

The  basic  solution  (i.e.,  for  cd  =  0)  is,  of  course, 


(21b) 


(22) 
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Since  these  are  constants ,  this  problem  is  similar  to  that  of  ordinary 
linearized  theory  for  steady  supersonic  flow  about  slender  bodies. 


The  introduction  of  a  and  (3  is  somewhat  arbitrary  as  long  as 
equations  (13)  are  satisfied;  these  parameters  may  be  conveniently 
taken  so  that  the  two  families  of  characteristics  of  the  basic  solution 
are  given  by 


,  x  =  aQt  +  p 

► 

x  =  -aQt  +  a 

so  that 

x  = 

2 

> 

4.  _  «  -  3 

"  ~^T 


(23) 


(24) 


For  method  B,  x  and  t  are  replaced  by  X  and  T  in  equations  (23) 
and  (24) . 


When  (a  has  a  finite  value,  the  xt-plane  may  be  divided  into 
regions  by  the  characteristics  which  pass  through  the  points  (xq,0) 
and  (Xp,0)  as  shown  in  figure  2.  The  solution  in  region  I  is 

unaffected  by  the  walls  at  the  ends.  In  region  II  the  effect  of  the 
wall  at  x  =  xq  only  is  felt;  whereas,  in  region  III  that  at  x  =  x^ 

only  is  felt.  Both,  walls  affect  the  flow  in  region  IV. 


The  solution  in  region  I  is  found  first.  (in  appendix  B  the 
solution  in  region  I  is  found  in  closed  form. )  With  method  A,  equa¬ 
tions  (22)  and  (24)  are  substituted  into  equations  (15).  The  first 
equation  of  equations  (15)  is  integrated  with  respect  to  a,  the  second 
with  respect  to  3,  so  that 


ui +  ± 


or  +  2ap  +  f(p) 


P) 


u,  -  Xa,  =  -  -i 
1  1  8a, 


0 


2ap  +  P^  +  g(a) 


(25) 
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where  f((3)  and  g(a)  are  arbitrary  functions  subject  to  the  boundary 
conditions.  These  functions  are  determined  by  noting  that  at  t  =  0, 
a  =  (3,  by  equations  (23) ,  and  that  u-j_  =  -  0  by  equations  (21a). 

It  is  found  that  f ( (3 )  =  -3(32  and  g(a)  =  -3a2.  By  substituting  these 
expressions  into  equations  (25)  and  solving  for  u-^  and  a]_,  the  fol¬ 
lowing  expressions  are  obtained: 


-1  -  xt 


ts 

\  2  )  21 


an  =  -  — — 
1  21a, 


l0 


(26) 


Equations  (l 6)  can  now  be  solved  in  the  same  manner.  Thus, 

^  /~  1  o\  /«  o\  3 


Up  = 


a0  = 


3a03 


-  -i*t3 


21  + 


2  =  2la2a03 


3  fa.  -  P\1*  <2X  +  3>a0  tlt 

=  ~ To2 


(27) 


Substituting  these  expressions  into  equations  (9)  gives  as  the  solution 
in  region  I 


2  1  4  1 

u  =  0)  xt  -  •=•  cjd  xtJ  +  .  . 


-§u  =  1  _  -L  CD2t2  +  2X  +  .3.  +  .  .  . 

a0  21  24l2 


(28) 


The  solution  by  method  B  starts  with  equations  (19a),  which  are 
the  same  as  equations  (15)  in  this  case.  The  solution,  then,  is 


1  (a  -  b\2  =  ^0  m2 

21an  \  2  /  21 


A1  _ 


(29a) 


# 
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Now  that  these  values  are  known,  equations  (19b)  may  be  solved.  The 
arbitrary  functions  are  determined  by  making  lines  in  this  solution 
pass  through  the  same  points  on  the  x-axis  of  the  xt-plane  as  the  cor¬ 
responding  lines  of  the  basic  solution)  that  is,  =  Tq  =  0  when 
a  =  3.  Thus, 


T,  =  L±. -l/fa  -  =  4  +  1  m3 

1  6Xa  A  2  )  "  6* 


(29b) 


Solving  equations  (20)  in  the  same  way  gives 


U0  =  1/a  ±  3  \/a 

2  2  A 


_  2X  +  1  jjij.3 


A2  ~  ~ 


-A 3) 4  =  .  gx  +  1  a  Ti 
,3  \  2  ;  .  oJj  2  a0T 


(30a) 


X2  = 


+  V)  fa  -  P\4  =  ^  +  4  4 

lA  2  A  2  /  24X 


To  -  (4  i  1)(5X  +  7)^a  -  _  (X  +  1)(^X  +  7) 


120X  a. 


(30b) 


Substitution  into  the  series  (17)  gives  as  the  solution  for  region  I 


u  =  od2xt  +  ±  1  ofaxT3  +  .  .  . 

6\ 


=  1  -  —  _  21  +  1  ^4^4 


(31a) 
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x  =  X  +  -  co2XT2  +  p-  cnhCT^  +  .  .  . 

2  241 

t  =  T  +  ±±±  0>2t3  +  (^l)(5X,t,7)  ^5  + 
6x  120  X2 


(31b) 


The  variables  X  and  T  may  be  interpreted  simply  as  parameters 
of  the  solution.  Actually,  they  define  points  of  intersection  of  char¬ 
acteristics  in  the  basic  solution  corresponding  to  characteristics  which 
intersect  at  x,t  in  this-  solution.  In  any  case,  X  and  T  may  be 
eliminated  from  equations  (31)}  the  procedure  results  in  equations  (28) 
again. 

The  only  difference  in  the  solution  for  region  II  is  in  the  deter¬ 
mination  of  the  arbitrary  functions  resulting  from  integration.  Those 
functions  resulting  from  integration  with  respect  to  3  may  be  deter¬ 
mined  by  the  condition  that  the  solution  must  equal  that  for  region  I 
on  the  boundary  between  the  solutions ;  that  is,  when  3  =  Xq.  The  other 

functions  can  then  be  found  from  the  boundary  condition  given  by  equa¬ 
tions  (21b).  That  condition,  however,  does  not  apply  to  the  solution 
of  equations  (19b)  and  (20b).  For  these  equations,  the  boundary  con¬ 
dition  is  replaced  simply  by 


X  =  X  +  0)2Xjl  +  OD^Xg  +  '  *  •  =  Xq 


or 


.  =  0 


=  x0> 


Xi  «  x2 
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The  solution  found  by  method  A  is 


%(x  -  xo)  (a0t  +  xo)  -  ;^^x  -  xo)|8^o3t3  + 

3(31  -  7)xQaQ2t2  -  12^(l  +  4)xq  -  (l  +  2)xj  XQSgt 
|(17X.  +  31)xq  -  (51  +  7)x](x  -  x0)xQ.  +  .  .  . 


oo2  r 

=  1  "  2AT^°2t2  +  2X°aQt  "  2X0(X  "  xo]|  + 


(21  +  3)aQ^t1+  +  (31  +  5)x0a03t3  +  I2(x  -  x0)X()a02t2  + 


3^(1  +  3)x0  -  (51  +  7)xj  (x  -  x0)xoaot  - 
4(1  +  2)(4x0  -  x)(x  -  x0)2x0]  +  .  .  . 


14 


NACA  TN  2794 


and  by  method  B 


u  = 


CD 


XaQT  -  x0(a0T  -  X  +  x0) 


CD 


X, 


:0(a0T  -  X  +  x0)4 


U  24A.aQ3 


4(21  +  l)a02T2  -  (13X.  +  11)  (X  -  x0)a0T  - 


k(2X  +  l)Xan3T3  - 


(X  -  xq) 


_a  _  1  cd 
2Aa  ,r 


(X  -  1)X  +  (114  +  13)xr 


>  + 


>(33a) 


aQ2T2  +  2x0(a0T  -  X  +  x0) 


cd 


24A,2  aQ^ 


(2X  +  l)anV  + 


x0(a0T  -  X  +  xq)<{ 


(7A  +  5)a02T2  + 


(ill  +  13)xn  -  (54  +  T)X 


a0T  - 


2(X  -  x0) 


(54  +  7)xn  -  2(4  +  2)X 


>  +  •  • 


2  r- 


x  =  X  + 


CD 


2Aa, 


AXaQ2T2  -  x0(a0T  -  X  +  x0)(A,a0T  +  X  -  Xq) 


t  =  T  + 


(X  +  1  )cd‘ 
12A.a03 


.2  1- 


2a3T3  +  3xQ(a0T  -  X  +'x0)(a0T  +  X  -  x0) 


>  (33b) 


The  solutions  for  region  III  are  obtained  in  a  similar  way  and  are 
the  same  except  that  Xq  is  replaced  by  Xj  and  the  signs  of  some  terms 

are  reversed.  In  regions  II  and  III^  it  is  found  that  the  elimination 
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of  X  and  T  from  the  solution  by  method  B  results  again  in  the  solu¬ 
tion  by  method  A. 

Figure  3  shows  characteristic  diagrams  for  the  case  Xq  =  2xq, 

1  an 

=  o  v  >  an<^  ^  =  5  (l  =  1.4).  The  dashed  lines  are  the  character- 

2  x0 

istics  of  the  basic  solution  and,  therefore,  are  the  mathematical 
characteristics  used  in  the  higher-order  solutions;  these  lines  are  in 
the  xt-plane  for  method  A  and  in  the  XT-plane  for  method  B.  The  solid 
lines  are  the  true  characteristic  lines  of  the  problem  in  the  xt-plane. 
They  are  found  automatically  by  method  B.  Method  A,  as  described  so 
far,  does  not  give  the  true  characteristic  lines,  but  it  is  found  that 
the  solution  of  equations  (7),  after  the  substitution  of  the  values  of 
u  and  a  (for  example,  eqs.  (28)  for  region  I) ,  results  in  the  same 
characteristic  lines  as  have  been  found  by  method  B.  This  determination 
is  necessary  for  the  evaluation  of  the  flow  in  the  field,  as  illustrated 
in  figure  4  where  typical  distributions  of  velocity  and  speed  of  sound 
are  shown  for  a  given  time  and  in  figure  5  where  typical  variations  with 
time  at  a  given  value  of  x  are  shown.  In  figures  4(a)  and  5(a),  the 
results  found  by  method  A  without  the  determination  of  the  variation  of 
the  characteristics  are  shown.  That  is,  the  lines  PR  and  PS  of 
figure  3  are  used  as  the  boundaries  between  the  regions.  It  is  seen 
that  u  and  a  are  discontinuous  at  these  boundaries  (although  the 
quantity  u  -  Xa  is  continuous  along  PR  and  u  +  Xa  is  continuous 
along  PS).  In  figures  4(b)  and  5(b),  the  results  are  shown  when  the 
true  characteristics  of  the  problem,  the  lines  QR  and  QS  of  figure  3, 
are  taken  as  the  boundaries.  Method  B  automatically  gives  the  results 
of  figures  4(b)  and  5(h).  Parts  (a)  and  (b)  are  identical  except  where 
they  represent _ conditions  in  the  regions  between  PR  and  QR  and 
between  PS  and  QS  in  figure  3* 


DISCUSSION 


Some  aspects  of  the  preceding  example  Vrhich  are  due  to  the  fact 
that  U  and  A  are  constants  and  which  would  not  ordinarily  be  present 
should  be  noted,  (l)  The  unsolved  equations  of  method  A  are  unusually 
simple  in  comparison  with  those  of  method  B  since  the  terms  containing 
the  unknowns  explicitly  drop  out.  (2)  As  a  result,  the  solution  of 
the  equations  of'  the  first  approximation  is  the  same  by  both  methods. 

(3)  The  solution  in  method  B  for  Xq  and  Tq  affects  the  unknowns 

u  and  a  only  in  the  second  order  and  similarly  for  the  higher 
approximations . 

The  basic  solution  cannot  usually  be  found  analytically  but  must  be 
determined  by  a  step-by-step  procedure.  In  that  case,  the  approximations 
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must  also  be  done  step  by  step.  The  characteristic  equations  would  be 
written  in  a  finite  difference  form  with  Aa  and  A3  cancelled  out. 
Thus,  a  and  3  would  not  be  used  explicitly.  With  method  B,  the 
calculations  would  be  carried  out  in  the  XT -plane  where  the  character¬ 
istic  directions  remain  constant  as  they  do  in  the  xt-plane  in  method  A. 
The  transfer  to  the  xt-plane,  which  consists  of  merely  adding  the  terms 
of  equations  (ljb),  can  be  delayed  until  the  desired  number  of  approxi¬ 
mations  has  been  made. 

From  the  preceding  simple  example,  some  inferences  can  be  drawn 
concerning  the  application  of  the  two  methods  to  other  basic  equations. 
The  unknowns  do  not  occur  explicitly  (although  they  do  appear  as  deriv¬ 
atives)  in  equations  (19a)  and  (20a),  but,  if  they  occur  thus  in  the 
basic  characteristic  equations,  corresponding  to  equations  (18a),  they 
will  occur  in  the  equations  corresponding  to  (19&)  and  (20a).  If  the 
original  independent  variables  (X  and  T  in  the  given  example)  appear 
in  the  basic  characteristic  equations,  the  four  equations  in  each  approx¬ 
imation  of  method  B  must  all  be  solved  simultaneously,  rather  than  as 
two  pairs  of  equations  as  in  the  example.  This  complication  is  a  serious 
disadvantage  of  method  B  and  may  rule  it  out  completely  in  such  cases. 

The  relative  simplicity  of  the  equations  of  method  B,  taken  one  by 
one,  is  probably  quite  general.  The  appearance  of  derivatives  of  lower 
approximations  with  respect  to  the  original  independent  variables  in  the 
equations  of  method  A  is  to  be  expected  as  long  as  the  primary  equations 
are  nonlinear,  whereas  they  cannot  appear  in  method  B. 

If  an  application  is  made  to  the  two-dimensional  supersonic  flow 
about  a  single  body,  for  example,  and  conditions  on  the  surface  of  the 
body  only  are  desired,  the  calculation  of  the  true  characteristics  is 
superfluous.  In  cases  such  as  this,  method  A  has  the  advantage  of  having 
two  less  equations  to  solve  in  each  approximation.  However,  the  pre¬ 
ceding  problem  indicates  that  the  true  position  of  characteristics 
through  points  of  discontinuity  in  the  boundary  conditions  is  necessary 
for  the  determination  of  the  flow  in  the  field  and  for  the  determination 
of  conditions  on  another  boundary  if  such  characteristics  intersect  it. 

# 

This  necessity  for  correcting  characteristics  as  it  applies  to 
linearized  theory  for  steady  supersonic  flow  is  discussed  in  reference  6. 
There,  as  in  the  present  example,  what  is  here  called  the  basic  solution 
consists  of  constants  so  that  this  effect  is  felt  only  in  the  second  and 
higher  approximations.  In  more  general  cases  it  would  be  felt  in  the 
first  approximation. 

In  problems  where  the  perturbation  involves  a  change  in  the  position 
of  the  boundary,  method  B  has  the  advantage  that  the  position  is  constant 
in  the  plane  of  the  original  independent  variables  of  the  basic  equations 
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(for  example,  the  XT -plane  in  the  present  problem)  where  the  solution 
is  carried  out. 

When  the  perturbation  involves  the  addition  of  a  new  dependent 
variable,  it  can  be  handled  in  much  the  same  way  as  the  entropy  is 
handled  in  this  case  (appendix  A). 


CONCLUDING  REMARKS 


Two  methods  of  using  the  concept  of  linearized  characteristics  are 
derived  for  the  one -dimensional  unsteady  flow  in  a  rotated  tube.  A  com¬ 
parison  of  both  methods  has  shown  that  they  give  the  same  results  but 
there  are  significant  differences  in  their  application. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  July  11,  1952. 
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APPENDIX  A 

DEVELOPMENT  OF  EQUATIONS  FOR  VARIABLE  ENTROPY 


When  the  entropy  is  not  constant  throughout  the  fluid,  equations  (5) 
are  replaced  by 


Xa-fc  +  aux  +  ^uax  ""0 

u-t  +  uux  +  Xaax  -  y  a2sx  -  u>2x  =  0^ 


(34) 


where  s  is  the  entropy  increment  divided  by  the  gas  constant.  As  long 
as  the  flow  is  without  shocks,  the  entropy  obeys  the  relation 

s-t  +  usx  =  0  (35) 

The  characteristic  equations  become 


Ua  +  y  asa 

-  *-ap  +  y  as3 

xa  =  (u  +  a)  to, 

Xp  =  (u  -  a)tp  • 


-  o)2xta  =  0 

p 

-  a)  xtp  =  0 


(36) 


and  equation  (35)  may  be  replaced  by 

satp  +  spta  =0  (37) 

(See  ref.  3  for  a  derivation  of  the  preceding  equations  for  o>  =  0.) 
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Writing 


s(x,t)  =  S  +  a>2sp  +  a)^s2  +  .  .  . 


substituting  this  expression  along  with  equations  (9)  into  equations  (34) 
and.  ( 35) >  and  converting  to  characteristic  form  give  the  equations  of 
method  A  as 


Ua  +  lAa  -  —  AScc,  -  0  j 

7 


Up  -  “p  +  7  *sp  =.°  > 


Satp  +  S)3ta  “  0 


(38) 


“h 


la  + 

Xala  +  ( 

[vx  *  wx)to(ul  +  al) 

■IP  " 

Xapp  +  | 

't  -  "^(“l  -  al) 

+  ?[AS1P 

la'tp 

+  sipta 

+  2UlSxtatp  =  0 

J 

(39) 


u2a  +  la2a  +  (ux  +  7AX)  ta  (v^  +  a2^  +  (ulx  +  la^tQ,^  +  a^ 


"h 


yjAs2a  +  (2al  +  ul)sla  +  (2a2  +  +  (U1  +  al)2sxto|  "  0 

U2p  -  Xa2p  +  (ux  -  Mx)tp^u2  -  a2)  +  (ulx  -  ^apx)tp  ^  -  a-^  +  ^  ^ 

i[As2P  +  (2ap  -  u^spp  +  (2a2  -  u2)sp  -  (u±  -  a^2SxtpJ  =  0 

s2at3  +  s2^a  +  2(ulslx  +  '^xJV'p  =  0 
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and  so  forth.  Equations  ( 13)  still  give  the  mathematical  characteristic 
direction  for  each  approximation.  If  the  true  characteristics  are  to  be 
found,  equations  (7)  still  apply. 

In  order  to  obtain  the  equations  for  method  B,  the  expression 


s(a,0)  = 

along  with  equations  ( IT )  is 
so  that 


substituted  into  equations  (36)  and  ( 37 ) > 


Uq,  +  XAq,  -  —  AS(x  -  0 
Up  -  AAp  +  i  ASp  =.  0 

Xa  -  (U  +  A)Ta  =  0 
Xp  -  (U  -  A)Tp  =  0 
Sa^P  +  SpTa  =  0 


(41) 


Ula  +  ^la  "  j  ^l^a  +  ^la)  ~  ^A 
uip  -  XA]_p  +  y (AjSp  +  AS-|_p^  =  XTp 
Xia  -  (U  +  A)  Tla  =  (Ui  +  Ax)Ta 
xip  -  (U  -  A)T1P  =  (ux  -  AjTp 
slaTP  +  saTlp  +  SppTa  +  SpTia  = 
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U2a  +  ^2a  "  ^(^2^a  +  ^l^la  +  ^®2a)  +  -^la 

U2(3  -  XA2(3  +  y(A2Sp  +  A^S^  +  AS2pj  =  XxTp  +  XT]^ 

X2a  -  (U  +  A)T2a  =  (UX  +  AjT^  +  (U2  +  A2)Ta 

X2p  -  (U  -  A)T2p  =  (Ui  -  Ai)t1P  +  (U2  -  A2)Tp 

S2c£Tp  +  SlaTip  +  SaT2p  +  S2j3Ta  +  SipTla  +  SpT2a  =  0 


and.  so  forth. 

It  is  to  be  noted  that,  when  the  entropy  of  the  primary  solution  S 
is  a  constant,  the  equations  of  both  methods  simplify  considerably.  For 
each  approximation  in  this  case,  the  equations  for  the  entropy  may  be 
solved  first,  then  the  equations  for  the  velocity  and  speed  of  sound, 
and  then  (in  method  B)  the  characteristic  direction  equations.  Other¬ 
wise,  all  three  equations  (in  method  A)  or  all  five  equations  (in 
method  B)  in  each  approximation  must  be  solved  simultaneously. 
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